Priedel oscillations due to Fermi arcs in Weyl semimetals 
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Weyl semimetals harbor unusual surface states known as Fermi arcs, which are essentially dis- 
joint segments of a two dimensional Fermi surface. We describe a prescription for obtaining Fermi 
arcs of arbitrary shape and connectivity by stacking alternate two dimensional electron and hole 
Fermi surfaces and adding suitable interlayer coupling. Using this prescription, we compute the 
local density of states - a quantity directly relevant to scanning tunneling microscopy - on a Weyl 
semimetal surface in the presence of a point scatterer and present results for a particular model 
that is expected to apply to pyrochlore iridate Weyl semimetals. For thin samples, Fermi arcs on 
opposite surfaces conspire to allow nested backscattering, resulting in strong Friedel oscillations on 
the surface. These oscillations die out as the sample thickness is increased and Fermi arcs from 
the bottom surface retreat and weak oscillations, due to scattering between the top surface Fermi 
arcs alone, survive. The surface spectral function - accessible to photoemission experiments - is 
also computed. In the thermodynamic limit, this calculation can be done analytically and separate 
contributions from the Fermi arcs and the bulk states can be seen. 
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Weyl semimetals (WSMs) are rapidly gaining popular- 
ity [l|43| as a new, gapless topological phase of matter, 
as opposed to topological insulators, which are gapped. 
A WSM is defined as a phase that has a pair of non- 
degenerate bands touching at a certain number of points 
in its Brillouin zone. Each such point or Weyl node 
has a chirality or a handedness; very general conditions 
constrain the right- and the left-handed Weyl nodes to 
be equal in number 0, S|- Near the nodes, the Hamil- 
tonian resembles that of Weyl fermions well-known in 
high-energy physics. These nodes are topologically sta- 
ble as long as translational symmetry is conserved, and 
can only be destroyed by annihilating them in pairs. Sev- 
eral theoretical proposals for realizing WSMs now exist 
in the literature 

@, S WSMs have already been 

predicted to exhibit several interesting bulk properties, 
ranging from unusual quantum hall effects fl3l [T3 | to 
various effects that rely on a 3D chiral anomaly present 
in this phase 
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-20J. Preliminary bulk transport 
studies of WSMs have been performed both theoretically 
0, HH, l22l|, as well as experimentally in some candidate 
materials I; 
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A remarkable feature of WSMs is the existence of un- 
conventional surface states known as Fermi arcs (FAs). 
These FAs are of a different origin from the FAs that exist 
in cuprate superconductors. A FA on a WSM is essen- 
tially a segment of a 2D Fermi surface (FS) that connects 
the projections of a pair of bulk Weyl nodes of opposite 
chiralities onto the surface Brillouin zone [2j. Although 
FAs always connect Weyl nodes of opposite chiralities, 
their exact shapes and connectivities depend strongly on 
the local boundary conditions. Such disconnected seg- 
ments of zero energy states cannot exist in isolated 2D 
systems, which must necessarily have a well-defined FS. 
A WSM in a slab geometry, however, is an isolated 2D 
system and indeed, FAs on opposite surfaces together do 
form a well-defined 2D FS. A natural question to there- 
fore ask is, "what signatures does this unusual Fermi sur- 



face have in scanning tunneling microscopy (STM) and 
angle-resolved photoemission spectroscopy (ARPES) - 
two common techniques that can probe surface states 
directly?" 

In this work, we answer this question by computing the 
local density of states (LDOS) on the surface of WSM, 
P sur f(f, E). in the presence of a point scatterer on the 
surface as well as the surface spectral function for a clean 
system, A ^(k, E). We apply our results to the iridates, 
A 2 Ir 2 7 , S l =Y, Eu, which are predicted to be WSMs 
di|. Both 

/'surf an< ^ ^surf evo ^ ve as ^ ne sam ple thick- 
ness is increased, and the evolution is explained in terms 
of the amplitude of the FAs on the far surface diminish- 
ing on the near surface. The calculation is done using a 
prescription that can give FAs of arbitrary shape and con- 
nectivity and simultaneously generate the corresponding 
Weyl nodes in the bulk. The procedure, in a nutshell, 
entails stacking electron and hole FSs alternately, and 
gapping them out pairwise via interlayer couplings that 
are designed to leave the desired FAs on the end layers. 
The resulting Hamiltonian is of a simple tight-binding 
form, which allows us to calculate ^4g ur f analytically in 
the thermodynamic limit. This quantity can be directly 
measured by ARPES. 

FAs appear in two qualitatively distinct ways: (a) ei- 
ther FAs on opposite surfaces overlap, resulting in a gap- 
less semimetal, or, (b) FAs on opposite surface do not 
overlap, resulting in a 2D metal with a FS. The 2D parti- 
cle density in this metal is proportional to the FS area ac- 
cording to Luttinger's theorem, and lives predominantly 
on the surface. In general, however, some particles will 
leak into the bulk, but the bulk filling will typically be 
0(1/ L), where L is the slab thickness, and will thus van- 
ish in the thermodynamic limit: L — > 00. In the model 
presented here, (a) results when equal numbers of elec- 
tron and hole FSs are stacked while (b) is obtained when 
the total number of 2D FSs is odd. This is consistent 
with the statement made earlier that the FA structure 
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depends strongly on the boundary conditions, since peel- 
ing off a single layer interchanges (a) and (b). 

FAs, in principle, can be generated by: (i) starting with 
a bulk model with the desired number of Weyl nodes, (ii) 
discretizing it in real space in the finite direction, and (iii) 
applying suitable boundary conditions to obtain FAs of 
the desired structure. While this approach works in prin- 
ciple, it has several associated complications. For exam- 
ple, determining the boundary conditions that result in 
the desired connectivity of the FAs is non-trivial. For in- 
stance, a WSM with four Weyl nodes W* 2 a t momenta 
xQi 2- where x — i denotes the chirality of the Weyl 
node, has two pairs of FAs on any surface on which the 
projections of the Weyl points are distinct. These FAs 
can pair up the Weyl points in two qualitatively different 
ways: as (W^W^) and (W 2 + W 2 ~) on each surface, which 
is an (a)-type connectivity, or as (IFy W-f) and (W^W-T") 
on the top surface and as (W^W^) and {W^W 2 ) on the 
bottom surface, which falls in the (b) category. However, 
there is currently no general prescription for determin- 
ing the boundary conditions that give one or the other 
connectivity. Moreover, to our knowledge there is also 
no general prescription for deriving lattice models with 
arbitrary numbers and locations of Weyl nodes. These 
gaps in working methods are filled by our top-down ap- 
proach for generating FAs directly. Our approach should 
be useful to model FAs in real systems, where surface 
effects can bend the FAs and change their connectivity 
unpredictably. 

Layering prescription: We describe the prescrip- 
tion by considering the simplest WSM, which has just 
two Weyl nodes, at (k x , k y , k z ) = (fe, k z ) = (Ki,2, 0) and 
the FAs connect K.\ and K 2 along a segment S (S') on 
the z = \ {z — L) surface in the surface Brillouin zone, 
as shown in Fig [T] S — S' and S ^ S' correspond to the 
two qualitatively different situations (a) and (b) men- 
tioned earlier, and will be obtained by distinct boundary 
conditions. Generalization to more Weyl points and Weyl 
points away from the k z = plane is straightforward. 

We claim that this WSM is generated by the following 
Bloch Hamiltonian: 

L 

H k = XXJ -1 )*^*.* (i) 

z=l 
L-l 

T Zjh hz,kipz+i,k + h.c. 



where even (odd) L generates S = S' (S ^ S'), £k is a 
phenomenological function that vanishes along a contour 
C given by 
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S = S' 



= SUS' S^S' 



(2) 



and the interlayer coupling h z j~ — — ifc(Afc) if z is even 
(odd). If S — S', C can be chosen arbitrarily as long as it 
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Figure 1: Above: Layering prescription for obtaining FAs of 
arbitrary shape. Dotted (dashed) ellipses represent electron 
(hole) FSs, and solid red segments are the residual FAs on 
adding interlayer hoppings t*. and A&,. The horizontal black 
dashed line on the topmost layer separates regions with Afe > 
tk and Afe < t^. An even (odd) number of total layers gives 
identical (non-identical) FAs on the two surfaces, as shown on 
the left (right). Below: ID systems at fixed fe £ C under the 
influence of A and t in the extreme cases where the smaller 
hopping vanishes, for even (left) and odd (right) L. Filled 
(empty) circles denote a state on an electron (a hole) Fermi 
surface in the limit of decoupled layers. The ellipses enclose 
the states which get mutually gapped out by the hoppings. 



contains the entire segment S. The functions and 
are real, non-negative phenomenological functions that 
satisfy: 

tJ >Ak k£S (3) 
\<A k k e C/S{= S' if C = S U S") 

^ dictates that tk = exactly at fe = K.\ t2 . The 
fc-dependence of t and A away from C is unimportant 
for our purposes, and will be assumed to be negligible 
henceforth. We now justify the above claim. 

Bulk: If the interlayer couplings tk — A& = 0, then 
EI = Hk describes a stack of alternate non-degenerate 
electron and hole FSs. When tf. and A fc are turned on, 
these FSs get gapped out in pairs in the bulk. Indeed, 
the bulk Hamiltonian is 

#fc,fc! k = £ kO- z + (Afc - tk cos k z )a x + tk sin k z a y (4) 

where z is the layering direction, and is gapped ev- 
erywhere except at the desired Weyl points: (k,k z ) = 
(■Ki,2, 0), due to ([3]). Allowing tk and Afc to be negative 
or complex simply moves the Weyl points off the k z — 
plane, but this does not affect the shape of the FAs. Near 
the gapless points, if bulk rea jj zes t ne W e yl Hamiltonian: 



H 



bulk 

[p- V^kA <r z 

p±v F (K l )a z ■ 



[P ■ Vfc (A Ki - t Ki )] v x + [t Kl Pz] cr y 
>\\Vi<i x + A p z a y (5) 
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where p± = p ■ e r (Ki) and p\\ = p ■ e t (Ki) are momenta 
perpendicular and parallel to C (e r (k) and e t (k) are 2D 
unit vectors normal and tangential to C) , vp is the Fermi 
velocity of the 2D FSs and v t = V fc (A K . - t K .)-e t (Ki). 
In going to the second line, the variation of tk and Afc 
perpendicular to C has been assumed to be negligible, 
since it does not affect the shape of the FAs. Vi has op- 
posite signs at K\ and K 2 , ensuring that the Weyl nodes 
have opposite chirality. -B^ulk * s obviously unaffected by 
the boundary conditions at z — 1 and z — L for large L. 

Surface: The surface, however, strongly depends on 
the boundary conditions; in particular, it is qualitatively 
different for odd and even L. If L is odd, at each k e C, 
a state remains unpaired and hence, gapless, at z = 1 
(z = L) whenever Afc < tk (Afc > tk)- The gapless 
states at z = 1 (z = L) thus, trace out S (S 1 ). On the 
other hand, when L is even, both ends of a the chain at 
fixed k carry a gapless state when Afc < tk and neither 
end has gapless states when Afc > tk- In this case, the 
gapless states on both surfaces of the slab trace out S. 

Viewed differently, the ID system at fixed k g C and 
|A| 7^ \t\ is an insulator in the CII symmetry class, which 
is known to have a Z topological classification 2^, 2(| . 
While |A| > \t\ gives a trivial phase, |A| < \t\ is topolog- 
ical^ non-trivial with a zero mode at each end protected 
by a chiral symmetry, if the ID lattice has a whole num- 
ber of unit cells. These end states are nothing but the FA 
states at that k when L is even. As k is varied along C, 
the ID system undergoes a topological phase transition 
at K\ and K%. For odd L. there is always a state at 
one end of the chain, as show in Fig [T] This prescription 
is similar in spirit to node-pairing picture of Ref [27| for 
chiral topological insulators in three dimensions. Note 
that it is necessary to start with two sets of FSs, since a 
single FS cannot be destroyed perturbatively. 

Symmetry analysis: WSMs can only exist in systems in 
which at least one symmetry out of time-reversal symme- 
try (T) and inversion symmetry (I) is broken; the pres- 
ence of both would make each band doubly degenerate 
and give Dirac semi-metals instead with four-component 
fermions near at each node. Moreover, I-symmetric, T- 
breaking (T-symmetric, I-breaking) WSMs have an odd 
(even) number of pairs of Weyl nodes. In the current pic- 
ture, the breaking of these symmetries can be understood 
as follows. Let us assume £k — £-&; if this weren't true, 
both symmetries would broken from the outset. In gen- 
eral, tk and Afc are unrelated to t_fc and A_fc, in which 
case both symmetries would again be broken. However, 
T is preserved if tk = t-k arid Afc = A_fc, which can only 
happen if the number of points on C at which tk = Afc, 
and thus the number of Weyl nodes, is an integer multiple 
of four. On the other hand, inversion about a particular 
layer interchanges t and A; thus, tk = A_fc preserves 
this inversion symmetry. This condition requires tk — Afc 
to change sign twice an odd number of times along C, 
giving an odd number of Weyl node pairs. 




(a) 



(b) 



(c) 



Figure 2: (a) FAs in the surface Brillouin zone on the 111 
surface of iridate WSMs. Solid (dashed) lines denote FAs on 
the near (far) surface. Filled (empty) circles denote projec- 
tions of Weyl nodes of positive (negative) chirality in all the 
figures above, (b and c) Surface LDOS in arbitrary units due 
to the six FAs near the L point in the presence of a point 
scatterer on the surface (brown cross) for a thin sample (b) 
and a thick sample (c). Insets show the numerically computed 
surface spectral function at E — for the clean system, with 
darker colors representing larger values. The computation is 
done for the model, described in the text, which generates the 
six FAs near the L-point but not the remaining eighteen FAs 
near the Brillouin zone edges. 



LDOS results: Having described the procedure for 
obtaining FAs from a 2D limit, we demonstrate its util- 
ity by calculating the surface spectral function for a clean 
system and the surface LDOS in the presence of a point 
surface scatterer within a model that should be relevant 
to the pyrochlore iridates A2L2O7, A=Y, Eu, which are 
purported WSMs with 24 Weyl nodes [g. The lattice 
in the WSM phase has inversion symmetry as well as a 
threefold rotation symmetry R3 about the cubic Ill-axis, 
and there are six Weyl nodes related by these symme- 
tries near each of the four L points in the FCC Brillouin 
zone. Additionally, the lattice also has a Dq symmetry, 
i.e., 7r/3 rotation about [111] followed by reflection in the 
perpendicular plane. This symmetry has an important 
implication for the FAs: if it is preserved in a slab ge- 
ometry, then the FAs will be as shown in Fig [2] (a). In 
particular, the six FAs near the center of the surface Bril- 
louin zone enclose an area, while the remaining eighteen 
overlap in pairs on the opposite surfaces. We note that 
Ref [3( also predicted Weyl semimetallic behavior in the 
above iridates, but with only 8 Weyl nodes. In this case, 
the hexagonal figure around the L point would collapse 
to a point. As we argue below, LDOS oscillations on the 
surface stem predominantly from the hexagonal figure; 
hence, the proposal of Ref if true, would imply no 
strong LDOS oscillations. 

We compute the surface LDOS for a model that has six 
FAs like the ones around the L point, as a function of the 
the sample thickness. The remaining eighteen FAs in the 
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iridates are expected to be destroyed by finite size effects 
for thin samples, while for thick samples backscattering 
occurs across the Brillouin zone and hence can give rise 
to LDOS oscillations only on the lattice scale. These 
oscillations are unlikely to be distinguishable from the 
electron density variations on this scale already present. 

The LDOS is calculated via the standard T-matrix 
formalism. Given the time-ordered Green's function for 
the clean system: Go(k,E) = (E — Hk) 1 and a scat- 
tering potential: U ZtZ >(x,y) — u5(x)S(y)S z ,i5z>,i, the T- 
matrix is given by T(lj) = (1 — UJ2k Go{k,E))~ U, in- 
dependent of momentum, since the scattering potential 
in momentum independent. Here, Go, U and T are all 
L x L matrices indexed by z. The full Green's func- 
tion in the presence of the impurity is G(fe, fe', E) = 
5 k . k >G {k,E) + G (k,E)T(u))G (k',E), and the LDOS 
on the z — 1 surface is related to the (1,1) element of its 
retarded cousin: p f(r,E) = — ilmG n (r, r, E + id), 

where G(r,r',E) = f k fc , e^ kr - h '- r '^G(k, k', E). For 

the calculation, we use £(k) = \Jk 2 + 2fc 6 cos 2 39 k — 1 
to generate the hexagonal figure and tk = to = 0.5, 
Afe = <o(l — cos36*fc) to obtain the FAs. Here, (k,0k) 
the polar coordinates of k. 

The results are presented in Fig|2](b) and (c) for E = 0. 
For thin samples, clear LDOS oscillations are seen in the 
horizontal direction as well as along the two equivalent 
directions related by 7r/3 rotation. The origin of these os- 
cillations becomes clear if one looks at A f (k,E = 0), 
displayed inset. Since the sample is thin, FA wavefunc- 



tions from the far surface have significant amplitude on 
the near surface, which allows backscattering to occur di- 
ametrically across the hexagon. The dominant backscat- 
tering processes are the ones involving the midpoints of 
the FAs, since the Fermi surface is nested here. On the 
other hand, as the sample thickness is increased, three of 
the six FAs retreat to the far surface and backscattering is 
exponentially suppressed. The result is small variations 
in the LDOS arising from scattering between FAs solely 
on the top surface. Thus, the STM map has a distinct 
evolution with sample thickness which is characteristic of 
the FA structure in the iridate WSMs. 

Surface spectral function in clean thermody- 
namic limit: ^(k, E) was computed numerically in 
order to generate the insets of Fig O In the thermody- 
namic limit, however, this calculation can be done ana- 
lytically. Denoting the (1,1) element of the clean Green's 
function for a L-layer slab by GJk^ , it is straightforward 
to show, by explicitly evaluating the (1,1) cofactor of E— 

H k and using det = det(A)det(£> - CA~ l B), 

that 



(E-£ k )GH L) (k,E) 



E 2 - el - K~ tl (e - £ k )Gl\ L _ 2) (fc, E) 



(6) 



In the thermodynamic limit: L — > oo, GjL^fc, E) 



°"0(L-2) ! 



k,E) = g{k,E). Thus, 
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g(k,E) 



K( E -e k ) 



(E 2 -£ 2 k + t\ 



£ 2 _i_ j.2 
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Ul(E* 



(7) 
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The physical condition A^^(k,E) = 
— -Im [g(k, E + iS)] > fixes the sign in front of 
the square root. Clearly, 



A° smi (k,E) = S(E-£ k )^ 
1 



t 2 - A 2 



\t 2 - A 2 I 

\ z k ^k\ 



2t U E - £ k ) 



W (E 2 -£ 2 +t 2 k - A 2 )' - At 2 (E 2 E 2 k ) 



(8) 



The first line is non-zero only when t\ > A 2 , and has a 
sharp peak at E = Clearly, this represents the con- 
tribution to ^4g ur f from the FA. Whereas, the second line 

is non- vanishing when \E\ > \£k \ and tk — \/E 2 — £ | < 

tk + VE 2 - £', 



< 



These inequalities are satisfied 
in the region near the projection of the Weyl points onto 
the surface. Moreover, this contribution to ^4g ur f has 



no delta-function peak. Thus, it represents contributions 
form the bulk states near the Weyl nodes. The quantity 
^surf can ^ e directly measured by ARPES experiments. 

In conclusion, we have studied impurity-induced 
Friedel oscillations due to FAs in WSMs, focusing on the 
FA structure of the purported iridate WSMs, and ob- 
served their dependence on sample thickness. For thin 
samples, FAs on both surfaces collude to allow nested 
backscattering and hence produce strong LDOS oscilla- 
tions, whereas for thick samples, the FAs on the far sur- 
face do not reach the near surface and such backscatter- 
ing and the consequent LDOS oscillations are suppressed. 
The calculation is done by building the desired WSM 
and FA structure by stacking electron and hole Fermi 
surfaces and adding suitable interlayer hopping. Within 
this prescription, the surface spectral function for a clean 
system can be calculated analytically in the thermody- 
namic limit. 
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